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Abstract 

A simple analytic expression for the non-singlet structure function fj\fs 
is given. The expression is derived from the result of Ref. [1] obtained by 
low X resummation of the quark ladder diagrams in the double logarithmic 
approximation of perturbative QCD. 
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The small x behavior of the non-singlet structure functions in DIS are 

of crucial importance for the determination of the quark densities. As an 
example, a verification of the Gottfried sum rules requires an extrapolation 
of the experimental data to region of low x [2] which can be done using 
non-singlet structure functions. 

The small x behavior of the non-singlet structure functions /^vs was 
studied in Rcf. [1]. Using the method of infrared evolution equation 
(IREE) the pQCD double logarithmic contribution to f^s was calculated. 
The result was shown to differ dramatically from the one predicted by 
the low X approximation of the conventional DGLAP equation. The non- 
singlet structure function f^s found in Ref. [1] was presented in the form 
of a complex lo integral, the inverse Mellin transform of a given partial 
wave. Though one can use this result to make general estimates and read 
off the leading asymptotics, it is not really "user friendly" and hard for 
numerical implementations. 

Recently a detailed study of the quark ladder double logarithmic resum- 
mations was carried out in Ref. [3]. Several methods were analyzed and 
compared. In particular, a comparison of a direct diagrammatic resum- 
mation using Bcthc-Salpetcr equation with IREE leads to a new insight in 
the structure of the latter. A particular outcome of this comparison is a 
clear understanding of the way the lo integrations should be performed. 

In the present paper we apply the results of analysis of Rcf. [3] to the 
results of Ref. [1] and derive a closed form expression for /ns- Fortunately, 
as observed in Ref. [3], a unique choice of a closed lo integration path 
allows the integration to be performed giving rise to a resulting expression 
in terms of modified Bessel functions Ii^. This expression is most easy to 
handle in any further applications. 

The structure function f^s is proportional to the imaginary part of the 
forward photon scattering amplitude projected onto the flavor non-singlet 
state. In the double logarithmic approximation of pQCD this amounts 
for summing up of diagrams built of two t-channel quarks and s-channel 
gluons forming a ladder. The quark ladders can be resummed using the 
IREE in partial wave representation. 

We define the partial wave through the ansatz: 

M{Q', s) = J^.[-^y F{^, Q'/f^') . (0.2) 

The function F{lo, Q"^ / /j,"^) is a. partial wave representation of the s-channel 
contribution to the quark-photon elastic amplitude M. The imaginary 
part is obtained from the 7i-crossed amplitude which at high energies ap- 
proximately equals M{Q'^, —s). The scale ji is introduced as an auxiliary 
infrared cutoff parameter, the minimal quark transverse momentum. For 
a more detailed description of the method of IREE we refer the reader to 
the original work [4] as well as to some applications in Ref. [1, 5, 6]. 

The IREE for F is obtained by differentiating (0.2) with respect to 
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In/x^. The equation reads 

dF 1 

+ = ^/oH^(^,y) (0.3) 

with y = ln(Q^//i^), The function /o was introduced in Ref. [4] 

foiu) = Att^ {u - ^a;2 - co^) ; = 2a,Cf /tt (0.4) 

The function F{uj, y) solving the IREE (0.3) has the form [4, 1]: 

F(a;,y) = ^/o(a;)e-(--^°/8.')^ for y>0, (0.5) 

The coefficient Cq contains information about coupling of the photon to the 
ladder, Co = — 4:Tr ae„ie{A) ■ e{A'), with e being the photon polarization 
vectors. In the Born approximation the above choice of Cq corresponds to 
a scattering off a free quark with initial condition proportional to (5(1 — x). 
The s-channel part of the elastic amplitude M reads 

M(Q2, s, jl) = ^ I ^ (^-^y fo{co)exp[yfo{u;)/87r']. (0.6) 

The non-singlet structure function f^s is related to M via 

-ttCo/ats = elimM ~ -vre^^^ (0.7) 

and this results in the following expression 

= ^ y ^ (qi) u^fo{u^)exp[yfo{cv)/87r']. (0.8) 

The results above are essentially copied from Ref. [1]. What is our new 
observation is that we can proceed further in performing the u integrations 
in Eqs. (0.6) and (0.8). The procedure is very similar to the one described 
in the Appendix B of Ref. [3] . 

As was emphasized in Ref. [3] the w-integration path Cent goes around 
the square root branch cut from —loq to uq. 



LO^ — LOn 



(0.9) 

Denote by ^ = ln{s /Q^) and rj = Iji{s/ /j,'^). Let us introduce a new 
complex variable 

z = {u - ^a;2 _ uj^)/ujo. (0.10) 

The uj integral then turns into a contour integral in the z-plane, and the 
integration path encircles the origin z = 0. We obtain: 

dz (\ 



0+2-^1 \z 



— Z I g'?'^o2/2 + Ca;o/2z 



(0.11) 
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Co (^Io('^o\/^) - ^hii^oVlv) 

The expression (0.11) is very similar to the one found in Ref. [3] for the 
non-singlet structure function of photon. 
Finally, the expression for /ns is given by 



Ins = 



4 ^f^ 



dz 



2] nLJoz/2 + ^uio/2z 



(0.12) 



9^ Cf [Jo+ 27ri V^:^ 

^(^/i(-ov/e^)-^/3(a;oye^)). 

Eq. (0.12) is the main result of this paper. In a somewhat different context 
a very similar expression can be also found in Ref. [7]. The expression 
(0.12) is correct for positive ^ and r/. At x = 1 (^ = 0), the u integral is 
divergent, which corresponds to S{1 — x) initial condition. 

An additional comment is in order. Strictly speaking the result ob- 
tained above should be called a structure function of a free quark rather 
than of a proton. This is because a particular choice of the initial con- 
ditions was implemented. We believe, however, that in the Regge limit 
considered above factorization holds implying the experssion (0.12) is a 
good approximation for a proton structure modulo an overall normaliza- 
tion constant. An alternative approach to the problem based on the NLO 
DGLAP equation and a low scale convolution can be found in Ref. [8]. 

In the Regge limit Ins/Q^ = Inl/x » InQ'^/jj,'^ 3> 1 the leading 
asymptotics of the Bessel functions coincide and cancel. 

fNS- (IT-J--(^T^' (0.13) 



Up to the pre-exponential factor this asymptotics was correctly found in 
Ref. [1]. Exploring another limit InQ'^/fi^ ^ Inl/x S> 1 we reproduce 
the low X approximation of the DGLAP equation 



fj,S ^gU^oVMl/x) InQV^^^ (0.14) 

ln(l/a;) 

Note that only the first term (7i) contrubutes to this asymptotics. Eq. 
(0.12) gives an analytic expression for the non-singlet structure function 
/jV5 and it provides a smooth interpolation between two high energy limits. 
It is worth a comment, however, that the low x approximation of the 
DGLAP equation is reproduced with a constant coupling constant only. 
An extensive work on inclusion of the running coupling effects has been 
done in Ref. [9]. 
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